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Abstract 

The coupling of a neutral hybrid {1, 3, 5 . . .} ^ exotic particle (or current) to two 
neutral (hybrid) meson particles with the same J PC and J = is proved to be sub- 
leading to the usual large- N c QCD counting. The coupling of the same exotic particle 
to certain two - (hybrid) meson currents with the same J PC and J = is also sub- 
leading. The decay of a {1,3,5...} h hybrid particle to rjir , r/'ir , rfrj, 7y(1295)7r°, 
vr(1300) 7r°, 7](U40)tt° , a (980)° a or /o(980)<7 is sub-leading, assuming that these final 
state particles are (hybrid) mesons in the limit of large iV c . 
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1 Introduction 

States of Quantum Chromodynamics (QCD) can definitively be said not to be conventional 
mesons when these states have exotic J PC , which cannot be constructed for conventional 
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mesons in the quark model, or equivalently, cannot be built from local currents with only a 
quark and an antiquark field. Here J denotes the internal angular momentum, P (parity) 
the reflection through the origin and C (charge conjugation) particle-antiparticle exchange. 
These are conserved quantum numbers of QCD. 

With the experimental discovery of isovector J PC exotics, the question of their interpretation 
has come into focus. QCD with a large number of colours N c offers a systematic expansion 
in 1/N C with considerable phenomenological success [1, 2], which can address this question. 
This is because a glueball (built from only gluons) and a (hybrid) meson (quark-antiquark 
with additional gluons) do not mix in large- iV c [1]. Furthermore, four-quark states (two 
quark-antiquark pairs) are absent [2]. In large- N c the isovector J exotics must therefore 
be hybrid mesons, as glueballs are isoscalar. Here it is proved for the first time that certain 
decays of hybrid mesons that are allowed by the conserved quantum numbers of QCD are 
sub-leading to their usual large- N c counting, providing a consistency check for the hybrid 
nature of the state. 

Selection rules for J PC exotic hybrid decay amplitudes, e.g. the amplitude for J PC = 1~ + 
hybrid particle — > r]7r, 77V, were noticed in non-field theoretic analyses [3]. In QCD it was 
found that these selection rules are really properties of certain three-point Green's func- 
tions [4, 5]. The first attempt to obtain hadronic properties from the Green's functions [4] 
contained some errors [5]. These properties were subsequently extracted in finite- N c QCD, 
e.g. the physical N c = 3 [5]. The properties were of limited physical relevance since they 
pertained to the coupling of currents to particles, e.g. a 1 _+ hybrid current to r/n. Also, 
for technical reasons, the scope of the deductions was limited. As will be seen below, these 
reasons disappear in large-iV c , because three and more particles do not contribute. The 
large- N c treatment of the results of Ref. [5] is the subject of this Paper. The results of this 
Paper can be comprehended by only reading Section 4, which also explicates the experimen- 
tal consequences of this Paper. Section 2 proves two results for the coupling of particles to 
currents by first proving two preliminary results. Section 3 uses the results of the previous 
section to prove a result for the coupling of particles to particles. 



2 Coupling of currents to particles 



The decay of a J PC = {1, 3, 5 . . .}~ + particle to two identical J = particles vanishes by 
Bose symmetry, because the final state particles are in an odd partial wave. The analogous 
statement for a Green's function built from a J PC = {1, 3, 5 . . .}~ + current and two J = 
currents is that the "identical current" part, or symmetric part, of the Green's function van- 
ishes [5]. The OZI rule allowed contributions to the Green's function only has a symmetric 
part [5], so that they do not contribute to the Green's function. The expression that does 
not contain an OZI rule allowed contribution is 



f°° dte iEt 6 v f d 3 x dV i(p -" p ' y) (0\B(x,t)C(y,t)A^0)\0) 

J —oo J 

= J2 (2?r) 4 5 3 (Pn) S(E n - E) 6 p (0| ( / d 3 x e* px B(x, 0)) C(0) \n) (nK(0)|0> . (1) 
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This is proved in Eqs. 2, 3 and 14 of Ref. [5] with no approximations. The left-hand side 
(L.H.S.) of the equation contains the time integral and spatial Fourier transform of a three- 
point Green's function which describes the "decay" of A into B and C. The expression is 
in Minkowski (physical) space with E and p real numbers. The L.H.S. is expanded on the 
right-hand side (R.H.S.) by inserting an infinite set of asymptotic stable states n with energy 
E n and momentum p n in order to extract physical predictions. The delta functions indicate 
that the asymptotic states are at rest and have energy E. The gauge-invariant local currents 
B, C and have the flavour structure of a neutral (hybrid) meson (linear combinations 
of uu,dd, . . . quark fields), and can contain gluon fields [6]. The currents B and C both 
have the same colour-Dirac-derivative-gluon structure for a given flavour, a finite number 
of derivatives (when expanded as a power series) and J = 0. Also, the currents B(0) and 
C(0) have equal P and C. The current A M (0) is assumed to have P = — and odd J (with 
Lorentz indices denoted by //). Conservation of charge conjugation then implies that this 
current is J PC = {1, 3, 5 . . .}~ + exotic, so that it should contain at least one gluon field: a 
hybrid meson current. 

Eq. 1 would be of limited interest were it not for the fact that the action of the operator O p 
(containing a finite number of derivatives in powers of p) allowed the demonstration that 
the L.H.S. contains only OZI rule forbidden contributions, and hence is 0(1) to leading 
order in the large- N c power counting, as opposed to the usual 0(N C ). The remainder of this 



Paper exploits this behaviour of the L.H.S. and deduces the consequences for the R.H.S. 
The strategy is to keep only the leading contributions to the R.H.S. in large- N c , and then to 
equate to the L.H.S. It is shown in the remainder of this section that the leading contribution 
to the R.H.S. comes from either J PC = {1, 3, 5 . . .}~ + one-hybrid-meson states, or from two- 
hybrid) meson states with J = and the same J PC (Eq. 2). The R.H.S. contains a product 
of two matrix elements for each of the leading contributions. One of the matrix elements 
will be shown to have the usual large- N c counting. The other matrix element will be shown 
to have one order in N c lower counting than usual (Eqs. 3-4), in order that the L.H.S. has 
an order lower counting than usual, as required. 

On the R.H.S. the usual large- N c counting for (0|f?(x, 0) C(0)\n) is order \fW C) N c , y/N~ c or 1 
for n respectively a one -, two -, three - or four - (hybrid) meson asymptotic state [1, 7]. The 
counting for (n|A M (0)|0) is respectively \fN~ c , 1, or l/N c [1, 7]. The product of the 

countings of the two matrix elements is N c , N c , 1 and 1/N C respectively. If the asymptotic 
states contained glueballs the counting of the product will be lower than N c . Hence only 
one- and two-particle (hybrid) meson states contribute in large- N c , and they contribute at 
0(N C ), as they should to equal the usual counting of the L.H.S. Also, the one-particle states 
that contribute to (n|A M (0)|0) at 0{y/W c ) are only neutral hybrid mesons with the same 
J PC as the current ^4^(0) [7]. These states cannot be mesons because they are J PC exotic. 
The two-particle states that contribute to (0\B(pc, 0) C(0)\n) at 0(N C ) are only two neutral 
(hybrid) mesons [2] with the same J PC as the currents B(0) and C(0). It follows that only 
one-hybrid-meson and two - (hybrid) meson states contribute on the R.H.S. to leading order 
in large-iV c . Using this the R.H.S. of Eq. 1 can be simplified to read (Appendix A.l) 

2vr £ 5(m„ - E) O p <0| ( / d'x B(x, 0)) C(0) \a 0) (a 0\A,(0) |0> + -±- £ (1 - %) 
xK(E) f dn kl 6 P (0| ( / d 3 x e ip x 5(x,0)) C(0) l^h^h) (^h^k^A^O) 

J J k 1 +k 2 = {0,E) 

(2) 

where the first sum is over the one-hybrid-meson states a (implicitly including the different 
polarizations); and the second over the two - (hybrid) meson states o~\ and <72, in such a 
way that a particular two-particle state is summed over only in the permutation o\0~i (and 
not <7 2 (7i ) in order to avoid double counting. Throughout the text the convention is that 
non-boldfaced variables starting with k (and p, q, x, y, z) indicate four-vectors. 

The kinematical variable dependence of the one-particle terms in Eq. 2 is only on E and 



p. If the hybrid particles have a discrete spectrum, there would only be contributions for 
discrete values of E when E = m a . Because of the constraint k\ + k 2 = (0,E), the two- 
particle terms also only depend on the kinematical variables E and p. Because K(E) only 
has support above the two-particle threshold, the two-particle terms all vanish below the 
lowest threshold, and contains a (continuous in E) contribution from each two-particle state 
above its threshold. 

In each of the terms in Eq. 2, one of the matrix elements has the usual large- iV c counting. 
For the one-particle terms this is (cr 0)^4^(0) |0> . The reason is that, barring accidental 
cancellations, this matrix element has to have exactly its usual large- N c counting (0(y/Nc)), 
in order not to violate the counting for two-point Green's functions (the first preliminary 
proved in Appendix A. 2). For the two-particle terms in Eq. 2, (0| £>(x, 0)C(0) \a1k1a2k2) has 
the usual large- iV c counting (0(N C )). The reason is that, barring accidental cancellations, it 
has to be exactly 0(N C ), in order not to violate the counting for four-point Green's functions 
(the second preliminary proved in Appendix A. 3). 

Eqs. 1-2 can schematically be written as L.H.S. = R.H.S. = J2 m a mbm + J2m c md m , where the 
sums are over a for the one-particle terms J2m a mb m and over <Ti<T2^ki for the two-particle 
terms J2 m c md m - (As shown in Appendix B, the integral over f2 kl can be written as a finite 
sum due to a partial wave expansion). The sum over a is finite, since the one-particle states 
must have mass equal to E. The sum over <ti<t 2 is finite because there is a finite number 
of two-particle thresholds below E. If the L.H.S. is subleading in the large-iV c counting, 
the R.H.S. is, but this does not imply that a m b m is subleading, nor that c n d n is, because 
it is possible that both a m b m and c m d m have the usual large-iV c counting, and there is a 
cancellation between the terms yielding a subleading R.H.S. The arguments in Appendix B 
have the consequence that a m b m and c m d m are each subleading, i.e. 0(1). This follows by 
evaluating the equation L.H.S. = J2 m a mb m + J2 m c md m multiple times for different currents, 
so that a matrix equation is obtained which is then inverted. This can only be done if 
the number of terms on the R.H.S. is finite, as it is here. If three- or more-particle terms 
contributed on the R.H.S., as would generally be the case for the physical N c = 3, it is 
not clear that the number of terms on the R.H.S. is finite, and that the inversion can be 
done. For this reason, the fact that three or more particles are subleading in the large- N c 
expansion, is important to make progress in the derivation. The remainder of this section 
is devoted to two results, which although exhaustively derived in Appendix B, are mostly 



evident at this point. 

The first preliminary ((er 1 (0) 1 0) has to be 0(y/~Nc)) means that b m is 0(^fW~ c ). Together 
with a m b m is 0(1), this implies that a m is 0(1/ \fN~c), which yields most of the first result 
of the Paper that the coupling of currents to a particle 

(0\B(x,t)C(y,t)\a0)=O(-±=), (3) 

where its usual counting is 0(^/W c ) . This holds for a neutral on-shell hybrid meson particle a 
at rest with J PC = {1, 3, 5 . . .}~ + . Also, B and C are neutral gauge-invariant local (hybrid) 
meson currents at space-time positions x and y at equal time with flavour structure a linear 
combination of uu, dd, . . ., the same colour-Dirac-derivative-gluon structure for a given 
flavour, a finite number of derivatives and J = 0. The currents B(0) and C(0) should have 
equal P and C. 

The second preliminary ((0| -B(x, 0) C(0) |<Ti/ci a 2 k 2 ) has to be 0(iV c )) yields most of the fact 
that c m is 0(N C ). Since c m d m is 0(1), this implies that d m is 0(1/ N c ), which yields most 
of the second result that the coupling of particles to a current 

(0-^0-^^)10) = O(i-), (4) 

where its usual counting is 0(1)- This holds for neutral on-shell (hybrid) meson particles 
o~\ and o" 2 with identical J PC and J = 0, and with arbitrary four-momenta k\ and hi- Also, 
t4 m (z) is a neutral gauge-invariant local hybrid meson J PC = {1,3,5...}~ + current with 
Lorentz indices at space-time position z with flavour structure a linear combination of 
uu, dd, . . .. 



3 Coupling of particles to particles 

In this section the dependence on the current A^(z) is removed from the matrix element in 
Eq. 4 to obtain a result (Eq. 6) that does not depend on the current, but on the physically 
relevant T-matrix. Because the matrix element in Eq. 4 has one order in N c lower counting 
than the usual, the T-matrix element (Eq. 6) will have an order lower counting than usual. 
It is shown in Appendix A. 4 that 



(aik 1 a 2 k 2 out \An(z)\Q) - (a 1 k 1 a 2 k 2 in \An(z)\Q) = ^ 



i2n 

7K 



[a (h + k^mlA^lO) 



x 8{Jm\ + kf + J ml + k| - m a ) 



JV C (crifci^^l T \a {k\ + fc 2 )) 



(5) 



when the restriction to the rest frame, ki + k 2 = 0, applies. This equation states that when 
the difference of the coupling of remote future "out" states and remote past "in" states to 
the current is considered, valuable information about the physically relevant T-matrix is 
obtained: the third result that the coupling of a particle to particles 



(a 1 k 1 a 2 k 2 \T\a0) = O(^ I ) , 



(6) 



where its usual counting is 0(l/\/W c ). This holds for neutral on-shell (hybrid) meson 
particles <J\ and a 2 with J = 0, identical J PC , and four- momenta k\ and k 2 in the rest 
frame ki + k 2 = 0; and for a neutral on-shell hybrid meson particle a at rest with J PC = 
{1,3,5...}-+. 

Even though the third result is proved in Appendix B using techniques analogous to those 
used to derive the first two results, its plausibility can be verified by using Eq. 6 in con- 
junction with the first preliminary ((a 0|A M (z)|0) has to be 0(y/W c )), to obtain that the 
R.H.S. of Eq. 5 is 0(1/N C ), consistent with the L.H.S. given by the second result (Eq. 4) as 
0(1/N C ). 



4 Remarks 

The three results of the Paper are Eqs. 3, 4 and 6, including the discussion under each 
equation. These results are theorems of large- N c QCD field theory with no approximations, 
and are valid within the generic large- N c framework [1, 2, 7]. 

The first result (Eq. 3) implies that certain four-quark currents are not good interpolators 
for hybrid meson particles. This may have implications for Euclidean space lattice QCD, 
even though the result was derived only in Minkowski space. A special case of the second 
result (Eq. 4) was previously derived [5] for an r/ir asymptotic state for certain quark masses 
within a certain kinematical range. 



The third result (Eq. 6) is of direct experimental relevance. For example, the decay am- 
plitudes (couplings of a particle to particles) of a {1,3,5...}~ + hybrid to r/n , rj'ir ,r]'r], 
r7(1295)7T°, 7r(1300) 7r°, ?7(1440)7r , a (980)V or / (980)<r are 0(1/nI), while the usual 
counting is 0(l/y/~Nc), assuming that these final state particles are (hybrid) mesons in the 
limit of large N c . Hence the widths of these decays are 1/iV^ suppressed with respect to 
their usual counting. This is the same suppression that large- N c predicts for decays forbid- 
den by the OZI rule [1], implying that the suppressions predicted here should be similar 
phenomenologically. The selection rule is most useful when OZI allowed decay is expected 
to be important in the absence of the selection rule. In the example above, this is true for 
a hybrid composed dominantly of uu and dd. An the other hand, it can be deduced from 

3 

Eq. 6 that the coupling of a 1~ + hybrid to r) c r) is 0(1/Nc), but this is less useful as the OZI 
allowed coupling (of the cc component of the hybrid to i] c (cc) and the cc component of the 
rj) is not expected to be important. Interestingly, even in the unlikely case where the rj' or 
a is a pure glueball in the limit of large N c , the decay amplitude of {1, 3, 5 . . .}~ + hybrids 
to r]'-K ,r]'ri, a (980)V or / (980)<7 would be 0(1/N C ) [1], which is still subdominant to the 
usual counting. In case a is a meson-meson state in the limit of large N c , the predictions 
mentioned do not apply. Beside the 0~ + and ++ particles mentioned, examples can also 
be given of + ~ and exotic particles in the final states. The large-iV c selection rules, 
in contrast to the selection rules discussed in Section 1, also apply when both final state 
mesons do not have the same radial excitation, e.g the ?7(1295)7r , 7r(1300)°7r° and ?7(1440)7r 
final states. Assuming isospin symmetry the results can also be extended to charged states 
by use of the Wigner-Eckart theorem, as will now be done. 

Consider the decay of a 1~ + isovector hybrid with isospin symmetry. Decay to r/ir, rj'ir, 
77(1295)7r, r/(1440)7r and a (980)o", which is ordinarily important, is suppressed. The exper- 
imental 7r 1 (1600) [8] is a 1 h exotic isovector resonance. It has not been seen in i]tt. A 1~ + 
enhancement at 1.6 GeV has prominently been seen in 7/7T [8], although the branching ratio 
is not dominant if the enhancement is resonant (-B(7r 1 (1600) — > /itt)/-B(7t 1 (1600) — > rj'ir) = 
3.80±0.78 [9]). If the enhancement is dominantly non-resonant, as has been advocated [10], 
the branching ratio is very small. The decay 7Ti(1600) — > r/(1295)7r is found to be small rel- 
ative to /i7T in an analysis of the r]n + n^n^ final state [9], although an earlier report stated 
that 7r!(1600) was seen in fiir and r/(1295)7r at a similar magnitude in K + K°-k~it~ [11]. If 
7Ti(1600) is found to have a large branching ratio to r/'n, that would be inconsistent with 
large- N c expectations which are otherwise consistent with its being a hybrid meson [12]. As 



discussed above, decay to r/'n is large- N c suppressed when rj' is either a meson or a glueball 
in the limit of large N c , although the suppression is less when rj' is a glueball. Hence a 
sizable r/'n branching ratio can arise through a large glueball component of an rj' meson [3], 
which violates the large-7V c prediction that meson-glueball mixing is suppressed. The re- 
cently discovered ^(2000) has not been seen in r/n, r/'n and ??(1295)7r [9], consistent with 
its being a hybrid meson. 
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A Appendix: Diverse results 
A.l Derivation of Eq. 2 

For the one-particle states, VJ n = VJ CT / d 3 p a / (27r) 3 , and the momentum p CT = due to the 



momentum 5-function, so that E a = ym^ + p CT = m a because the particles in an asymp- 
totic state are on-shell. For the two-particle states VJ n = ^2 ai(72 (l — 5 ai(T2 /2) J d 3 k± / (2n) 3 
J d 3 k 2 I (2"7r) 3 , where the factor (1 — S ai(72 /2) is 1/2 for the phase space of identical particles. 
Substituting p„ = ki + k 2 and E n = + m| + y / k|+~m| in the phase space integration 



/ w? I S <2w)1 s3(p " } m ~ - E) f(h ' *>> = w / m ' ■ 

(7) 

where 



fci+fc 2 =(0,E) 



K(E) = j™ k 2 d|k!| <f(0c? + m 2 + ^/k 2 + m\-E) 



y / , , 4 (E - (mi + m 2 )1mi - m 2 ) 2 ) y (-E 2 - (mi + m 2 ) 2 )(£ 2 - (m x - m 2 ) 2 ) (8) 
if E 1 > mi + m 2 ; and K{E) vanishes if E < mi + m 2 . 

A. 2 First preliminary 

The following two-point function is O(N ) in the Feynman-diagrammatic large- N c count- 
ing [7] 

(0\A^ Xl )Au(x2)\0) = ]T (Ol^^Oln) (n\A v (x 2 )\0) , (9) 



and only one-(hybrid)-mesons n contribute at leading order [7]. Hence there must be a 
non-empty set of states n for which (0|A M (xi)|n) (n\A 1/ (x 2 )\Q) of 0(N C ). If A, //, v, x\ and 
x 2 are changed the set of states for which this is true may change. As these variables are 
changed, a specific state n should regularly be part of the set, since there is nothing special 
about it. Hence for a specific particle a with four-momentum p a , it must be possible to 
choose A, /i = v and x\ = x 2 = z such that (0|-A /1 (xi)|n) (n\A u (x 2 )\0) = \ {ap a \A^{z) |0)| 2 is 
0(N C ). This implies that {ap^A^z)^) is 0{y/~N^), as promised. 

A. 3 Second preliminary 

The following four-point function is 0(N%) [2] 

(0\B(x 1 )C(x 2 )B(x 3 )C(x 4 ) |0> = J2(0\B(x 1 )C(x 2 )\n) (n\B(x 3 ) C(x 4 )\0) , (10) 

n 

and only two -(hybrid) meson states n contribute at leading order [2]. Similar to the ar- 
gument for the first preliminary, it must be possible to choose B, C, x\ = x 3 = (x, 0), 
x 2 = x 4 = and a specific state \crikia 2 k 2 ) , such that (0\B(xi) C{x 2 )\n) (n\B(x 3 ) C(x 4 )|0) = 
I <0 1 5(x,0) C(0) K&i^)! 2 is C{Nl). Whence the promised result. 

A. 4 Derivation of Eq. 5 

All derivations so far left unspecified whether the asymptotic states were "in" or "out" 
states. Consider the specific case of "out" states, and insert a complete set of "in" states: 



(<TiA;i 0-2^2 out \A^(z)\fy =J2 777T^ {o-\k\o 2 k 2 out | erg in) (aq in jO) 

a J ( 27r ) 

+ (1 — -^p") y y (tTifci^fe out |<7^i<7 2 g 2 in) (ct^i cr 2 g 2 in |A M (z)|0) . (11) 

Restrict to ki + k 2 = in this subsection. The usual large- N c counting for (crikio- 2 k 2 \n) is 
order \/y/N~ c , 1 (no scattering), or \/y/N~ c for n respectively a one -, two - or three - (hybrid) 
meson state [1, 7]. The counting for {n\A^(z)\ti) is respectively \fW c , 1 or 1/v^Vc [1, 7]. 
The product of the countings of the two matrix elements is 1, 1 and 1/N C respectively. If 
the asymptotic states contained glueballs the counting of the product will be lower than 1. 



Hence only one - and two - (hybrid) meson states contribute in large-iV c , as indicated in 
Eq. 11. As before, the one-particle states that contribute are only neutral hybrid mesons 
with the same J PC as the current A^(0). 

Connection with the S-matrix is now made by using (mout|nin) = (min|S'|nin) [13]. 
In the second part of Eq. 11 write the two-body scattering (aikia 2 k 2 out \o' x q\o' 2 q 2 in) = 
(p\k\0 2 k 2 in \o' x q\a' 2 q 2 in) where it was used that only no-scattering occurs at 0(1) [1]. The 
latter overlap is simply an overlap between free bosonic states in the same basis. It can be 
evaluated [13] and equals 

(2tt) 6 ( 5 3 (k! - qi ) 5 aia[ 5 3 (k 2 - q 2 ) 6 a2(7 , + 5 3 (k x - q 2 ) 6 aia , 5 3 (k 2 - Ql ) 5 a2< ) . (12) 

In the first part in Eq. 11 introduce the T- matrix, defined as the transition from an "in" to 
an "out" state 

{a±kia 2 k 2 out | erg in) = i(27r) 4 5 4 (ki + k 2 — q) (a±kia 2 k 2 in \T\aqm) , (13) 

using that no-scattering does not contribute, given that this is the overlap of a two- with a 
one-particle state; and employing the definition of the S-matrix in terms of the "reduced" 
T-matrix [13]. Dropping the "in" label on the R.H.S. of Eq. 13, because (min|T|nin) = 
(mout |T|nout) [13], and substituting Eqs. 12-13 in Eq. 11, yield Eq. 5. 

B Appendix: Why each term in Eq. 2 is subleading 

This Appendix starts by proving that the angular integration in Eq. 2 can be written as 
a finite sum. The first two results of the Paper are then derived. (This is done directly, 
without first showing that each of the terms in Eq. 2 is subleading as discussed in the main 
text. However, once the first two results are established, it follows that each of the terms is 
subleading). The third result is subsequently derived. 

It will be convenient for the derivation to write the integral over the solid angle in Eq. 2 
as a finite sum. This can be done by performing a partial wave expansion of the overlap 
(cri/cicr 2 /i;2|^t(0)|0) by explicitly considering its Lorentz structure. It is a function of ki and 
k 2 , or equivalently of k\ — k 2 and k x + k 2 . First consider the Lorentz scalars that can be 



built from these two variables: (ki — k 2 ) 2 , (ki + k 2 ) 2 and (k\ — k 2 ) ■ (k\ + k 2 ). It is easily 
shown that the on-shell conditions k\ = m\ and k\ = m 2 imply that the latter two variables 
can be expressed in terms of the first variable. Hence the only independent Lorentz scalar 
is (k\ — k 2 ) 2 . Second consider the case where has J = 1, i.e. is a vector. The overlap can 
be written as a linear combination of times a Lorentz scalar, and k 2fJj times a Lorentz 
scalar, since this is the most general structure transforming like a vector. Denote the two 
Lorentz scalars by ((Jicr 2 \A(0))i((ki — k 2 ) 2 ), with % — 1,2. Define Ct(ki,k 2 ) = k ilM . Then 
the overlap 



M^IA^IO) =Y^^(h,k 2 ) (<7 lf 7 2 |A(0)); ((h-k,) 2 ) . (14) 

i 

It is evident that the procedure can be performed for arbitrary J, and that appropriate 
can always be constructed, with the partial wave i ranging over a finite number of integers. 
Eq. 14 is the promised partial wave expansion of the overlap for general J. The functions 
depend purely on kinematical variables and all dynamical information is contained in 
the scalar functions (cri<j 2 \A(0))i {{k\ — k 2 ) 2 ), which depend kinematically only on {k\ — k 2 ) 2 . 
When Eq. 14 is substituted in the two-particle terms of Eq. 2 the integral over the solid 
angle f2 kl can be written as a finite sum over i, as promised, since (ki — k 2 ) 2 does not depend 
on the solid angle. 

Rewrite Eqs. 1-2 as 

W BC = E M BCa V a + J2 M BC(7ia2i V aia2i , (15) 

(7=1 <71<X22 —1 

where the explicit dependence on the currents B and C are indicated. The L.H.S. of Eq. 1 
(divided by N c ) is 

W BC = ^r [°° dt e iEt 6 P / d 3 x d 3 y <0| £(x,t) C(y,t) A^O) |0) . (16) 

The one-particle matrix elements of Eq. 2 (divided by iV c ) are given by 

M BCa = Nc' 1 " 2tt 5{m a - E) O p (0| ( J d 3 x e ip ' x B(x,0)) C(0) |<r0) , (17) 
V a = N~ a ~^ (a 0^(0)10) , (18) 

where a is a real number specified below. The two-particle matrix elements of Eq. 2 (divided 

by N c ) are 




k!+k 2 =(0,E) 



(20) 



(19) 



using Eq. 14. As discussed in the main text, each side of Eq. 15 depends on the kinematical 
variables p and E. For a specific choice of these variables, there are iV s > one-particle 
states a contributing, and Nu > two-particle states and partial waves a±a 2 i contributing. 
It will be useful to think of the group of labels o\Oii as a single label. The remainder of the 
discussion is only of interest if it is not the case that iV s = N u = 0. 

From the second preliminary ((0| 5(x, 0) C(0) \cr1k1a2k2) has to be 0(N C )) it follows that M 
in Eq. 19 is < 0(1). This obtains by noting that O p is independent of colour, as is shown 
at the end of this Appendix, and that and K are purely kinematical functions with no 
colour dependence. The possibility that there are accidental cancellations in the various 
integrations, which could make M < 0(1), is incorporated by indicating that M < 0(1). 
The possibility of cancellations will be taken into account in the derivations below. By 
choosing a appropriately, M in Eq. 17 is defined to be exactly 0(1). Hence both M and M 
are < 0(1). Evaluate Eq. 15 for N^ + Nji different currents B, C. Since V and V in Eqs. 18 
and 20 are independent of the choice of currents B and C, this amounts to constructing a 
matrix equation W = M V. The + iVn dimensional column vector W is built from the 
evaluations of W B c fc> r different values of B, C; the (iV s + N u ) x (iV s + N n ) dimensional 
matrix M contains M BC a and M BC a X aw and the iV s + iV n dimensional column vector V is 
built from V a and V aia2i . Because both M and M are < 0(1) it follows that each entry of the 
matrix M is also < 0(1). This means that, barring accidental cancellations, the determinant 
of M, det M, which is a sum of products of the entries of M, is exactly 0(1), i.e. is non-zero 
and finite in the large- N c limit. Note that even if some of the entries of the matrix M are 
< 0(1), it is still possible for det M to be exactly 0(1). If det M < 0(1) the derivations 
below are invalid. This possibility can be excluded by an appropriate choice of currents. 
Since det M = 0(1) the inverse of M exists, and V = M^ 1 W. Since M" 1 = adj M / det M, 
where adj M is the adjoint matrix of M, which is a sum of products of entries of M, it follows 
that the entries of M _1 < 0(1). Noting from the main text that the L.H.S. of Eq. 1 is 0(1) 
[strictly speaking, it is < 0(1), since it is only known that OZI allowed 0(N C ) contributions 



are not present [5], and that the highest order OZI forbidden contribution is 0(1) [1] if there 
are no accidental cancellations], so that each entry of the vector W < 0(1/N C ), it follows 
from V = IVU 1 W that each entry of the vector V is < 0(1/N C ). This implies that V and 

V are both < 0(1/N C ). 

It is instructive to study the kinematical variable dependence of V in Eq. 20, which can 
only be relevant to the discussion if E is above the two-particle threshold m\ + m 2 , since 
the two-particle term in Eq. 2 only has support above this threshold. This, together with 
the constraint k\ + k 2 — (0, E) in Eq. 20, and the on-shell character of the particles, can 
be shown to imply that (k\ — k 2 ) 2 = 2 [m\ + m|) — E 2 G (— oo, (mi — m 2 ) 2 ]. Using 

V < 0{1/N C ), and Eqs. 14 and 20, it follows that (a 1 k 1 a 2 k 2 \A fl {0)\0) < 0(l/N c ) with the 
constraint that (k\ — k 2 ) 2 G (— oo, (mi — m 2 ) 2 ]. However, it is possible to show by only 
considering the on-shell nature of the particles, that the same constraint holds. Hence the 
constraint adds no new information, and is dropped henceforth. Thus 

(a.haMA^lO) = e^+^i^ha.hlA.mO) < 0(1/N C ) , (21) 

using space-time translational invariance A^(z) = e tP ' z A^(0) e~ tP ' z , with P u the QCD four- 
momentum operator. It is evident that equality in Eq. 21 would have been attained were it 
not for the possibility of accidental cancellations. These cancellations can be eliminated by 
an appropriate choice of the currents. Whence the result in Eq. 4. 

The observation that V < 0(1/N C ), together with the first preliminary ((a 0|A^(0)|0) has to 
be 0(t/Nc)), implies from Eq. 18 that a > 1. Since M was defined to be (9(1) this implies 
from Eq. 17 that (0| (/ d 3 x e ip ' x 5(x,0)) C(0) \a0) is < 0(l/y/7Q, noting that m a is 
Oil) [7]. Inverting the Fourier transform, it follows that (0|£(x,0) C(0)|<rO) < 0(l/y/7Q. 
Using space-time translational invariance analogous to Eq. 21 

(0|5(x,0 C(y,t) \aO) = e- im ° t (0| 5(x - y, 0) C(0) \a0) < 0(1/ Jn~ c ) . (22) 

From the same arguments as those below Eq. 21, the result in Eq. 3 is deduced. 

It remains to prove that Eq. 6 can be deduced from Eq. 5. The proof is analogous to the 
proof already given in this Appendix, and the various steps are outlined. The notation of 
the vectors and matrices will be the same except that the label BC will be replaced by the 
label A. Call the L.H.S. of Eq. 5 Wa- The first and second terms in long brackets on the 



R.H.S. of Eq. 5 are called Ma<, and V a respectively. Eq. 5 is then of the form of Eq. 15 with 
N n = 0. From the first preliminary ((a 0\A^(z)\0) has to be 0(y/N^)), it follows that Ma<j 
is 0(1). Evaluate Eq. 15 for iV s different currents A. This again gives a matrix equation 
W = M V. Barring accidental cancellations, det M = 0(1), so that its inverse exists. From 
Eq. 4 W < 0{1/N C ), and together with M" 1 < 0(1), the equation V = M" 1 W implies 
that V < 0(1/N C ). Along the same lines as before this establishes the result in Eq. 6. 

It is lastly outlined why O p does not depend on colour. This is done by following the deriva- 
tion of O p in the Appendix of Ref. [5], employing the notations of that reference. Colour 
appears when the (anti) commutators are evaluated, e.g. as 5 ab in {^(x, t), V^(y, t)} = 
— 5 a6 7^ • <9 x 5 3 (x — y). The colour and Dirac indices in the commutators are then contracted 
with the remaining quark and gluon fields and subsumed in / ^(x , y , z) . The construc- 
tion of O p only depends on the number of derivatives acting on 5 3 (x — y) and not on 
f^(x,y,z), making it independent of colour, as promised. Let's give an example of how 
this observation is used above. Suppose O p = d/dp, then in Eqs. 17 and 19 it occurs as 
Op J d 3 x exp(ip ■ x) g(x) = ix. J d 3 x exp(ip • x) g(x). It is evident that O p does not affect the 
large- N c counting of the function g. 



